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This  note  was  prepared  as  a  result  of  a  suggestion 
froml-'r.  'ard  F.  j.avidson,  (consultant  to  the  Chairman,  I'D'^C, 
that  a  jnore  ezact  method  could  be  developc^d  for  aajustin,"^  a 
curve  to  creep  data  than  those  usually  employed  or  the  one 
sug  e&ted  by  McVetty  in  ^i  rjcent  issue  of  the  Transactions 
of  the  American  Society  of  Mechanici.l  Engin^^ers.   It  is 
hoped  th^.  t  the  method  outlined  here  will  be  useful  to  investi- 
gators and  engineers  concerned  with  the  problems  of  materials 
under  creep  conditions. 
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FITTING  A  HYI^ERBOLIC  Sllffi  CURVE  TO  CREE^=  DATA. 
By  M.  Shiffmr.n  and  B.  Friedman 

Various  more  or  less  empirlct.l  formulas  have  been  used 
to  connect  the  rate  of  creep  of  a  metal  (at  a  constant  tem- 
perature) with  the  stress  applied.   The  power  function 
r  =  a  s    where   r   Is  the  rate  of  creep,   s   is  the  applied 
stress,  and   a,  n   are  constants,  has  frequently  been  tried 
because  of  its  simplicity.   Hov/ever,  this  formula  often  ler.ds 
to  very  incorrect  conclusions,  especially  v,jhen  extrapolating 
from  experimentally  determined  points  to  smrller  stresses.   It 
is  more  appropriate  to  use  a  relation  which  has  some  theoreti- 
cal justification,  provided  that  this  cm  be  done  without  too 
much  labor  and,  on  the  other  hand,  pives  a  better  fit.   The 
"hyperbolic  sine  law" 

(1)  r  =  a  slnh  bs 

1      2 
has  been  suggested  on  theoretical  grounds  by  Prandtl,  Nadal 

and  Eyring,   and  the  question  is  to  fit  this  formula  to  a  set 
of  experimentally  determined  points  by  finding  suitable  a   and 
b. 

This  problem  is  discussed  in  a  pai)er  by  McVetty  in  the 
Transactions  of  the  American  Society  of  Mechanical  Engineers, 
October  1945.   There  a  graphical  method  is  exhibited  for  de- 
termining the  constants   a,  b   in  the  case  of  two  frivcn  points. 
V/hen  there  arc  more  thrj:i  two  given  points,  McVetty  determines 
a,   b   for  all  pairs  of  points  and  selects  the  average  of  these 
as  the  appropriate   a,  b.   This  procedure  turns  out  to  be 
reasonrbly  accurate  only  if  thu  points  happen  to  fit  an  expo- 
nential curve  very  closely.  V/hile  this  is  the  case  in  McVetty' s 
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example,  it  need  not  rjw?  ys  be  true.   Si-;e  illustrative  dr.tr 

on  pfge  4  ,  tr.ken  from  the  Creep  Rrtc  Book. 

In  this  memorandum  other  possible  procedures  for  deter- 
mining a,   b   are  developed  which  are  valid  in  rny  case  where 
the  hyperbolic  sine  law  is  apv;licable,  and  are  easily  handled 
numerical ly . 

1.  Largo  Stresses.  In  the  first  place,  if   s   Is  large 

(say  bs  >  2)   the  equation  (1)  can  be  approximrted  by 

_  a  bs 
^   2  e 

or    (2)  log      r   =   log  -g  +   b«iog   O'S. 

This    is   a    strain-ht  line   on   semi-logarithm  paper,    and   the    con- 
stants     a,    b      can   be   determined  by   standard  methods.      While 
the    assumption   of   large      s      is    not    alwrys    correct,       this   happens 
to  be    the    case    in  McVetty's   Illustration,    and   in   fig.    1  we   have 
plotted   the   data   on  semi -log  paper.      The   data   fits    a   straight 


line    closely;    measuring   the   slope   and   intercept  yields    the   vp.lue 

J. 
624 


of   a  and  b.   The  results  r  re   b  =  -^^^r  >    ''^  -   .0021.   Thus, 


we  have  quickly  obtained 

r  =   .0021   sinh  ■— 

which  fits  the  experimi^ntr.l  points  parti cul-.rly  closely.  Even 
in  this  case  I,  the  procedure  here  suggested  might  have  advan- 
tages   over   ricVctty's. 

II.    Arbitrary   Stresses.      In   the   gcner'^1    case    the    experi- 
mentrl   datr    plotted   on    seni-logr  ri  thjn  p^.per   will   not   fit   a 
strr.  ight    line    closely   but  will    desc^^nd  more    rapidly   for    smrller 

values   of      s.       Compare   v;ith   the   hyperbolic    sine;    curve   given   in 
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fig.  2.   A  least  squrro  criterion  C'n  bo  used  to  dutcrmino 
the  best  a,   b. 

Replace    r  =  a  sinh  bs 

(3)  log   r  =  A  +  log  sinh  bs   ,  whore  A  =  log  a, 
and  apply  the  least  squrrcs  criterion  to  (3),  that  is,  deter- 
mine A,   b   so  as  to  minimize 

(4)  .^    (log  r.  -  A  -  log  sinh  bs  l""   . 
i=l         ^  ^ 

This   seems   necessary   since    the    stresses      s.      are   deter- 
mined  r-ccur;itely   Vi/hile    the    experimental    error   of   the    creep 
rates      r.       are   undoubtedly   of    the   percentrge    type,    rather    thrn 
the   absolute    type.      This    justifies   the   use   of   log      r.      in   plrcc 

of   r.  . 

1 

The  equations  obtcined  by  minimizing  (4)  con t:  in  A  rnd  b 
1»  a  transcendental  mrnner.   But  this  set  of  equr.tions  cm  easi- 
ly be  solved  by  a  method  of  iteration,  of  v;hich  the  first  step 
will  usuf'lly  suffice. 

The  procedure  to  follow  in  mrking  the  approximation  will 
be  described;  an  example  v/ill  be  worked  out;  and  then  the  de- 
rivation of  the  oqup-tions  used  will  be  given  at  the  end: 

Pro  cedure  : 


First,  plot  the  experimental  points  on  a  semi-logarithm 
paper  rnd  determine  graphically  a  straight  line  which  approxi- 
mately fits  the  points.   The  slope  of  this  str.;ight  lino  will 
yield  a  first  approximation  b    to  the  vrlue  of   b.   The  meth- 
od is  the  sr.me  as  in  part  I. 
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The  first  approximation  /    to  the  value  of  A   is 
obtained  from 


(5)  x\q  =  -  .^^  (log  r^  -  log  slnh  h^s^). 

1—1 

The    corrections    Ab,    A  A      to  be   added   to   these   values 

b   ,      A        are   found  as  follows: 
o  o 

Define    the    "weights"   w.      by 

w.    =   s ,      coth     b   s, . 

11  o     1 

Let     w     be    the    average   of    these  weights,    w  =  -     ^   w..    Then 

y  (w.    -   w)    (log   r.    -   log   sinh  b    s.) 

(6)  Ab   =  : i ^^- 2~^-^ 

.4343      X     2.  (^"^1    -   ™") 

and 

(7)  A  A   =   -.4343  w.  Ab. 

The    quantity    ,4343   is    log      e. 
Example  : 

Consider    the    following   data    taken  from  the    Creep   Data 
Handbook,    1939,    page   245,    for   a  wrought    iron   bar    at   1200°  P 
under    stress    for   1,000   hours. 

s(lb./in.^) vi%  per   10,000   hrs.) 

400  .39 

750  .90 

1250  1.59 

1750  3.32 

2500  3.75 

Plotting    these   points    on    semi-log   paper    (fig.    2)    and 

selecting   a    suitable    straight   line,    b        can  be   determined   from 

the    slope   of  the   line.      'Ne    find   the    slope    to   be    .0006   or 


b        log      e   =    .0060  and     b^     =    .0014. 
o  o 
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The    remaining    calculations    can   be    set  up   as    ■f'ollows 


log   r 


400  .39  -.4089 

750  .90  -.0458 

1250  1.59  .2014 

1750  3.32  .5211 

2500  8.75  .9420 


-  -    2 

b   s  coth  b   s  s    coth  b   s        w-w  (w-w) 

o  0  o  ^  ' 


Sunr 


Average 


.56 

1.97 

788 

-684 

.468    X   10^ 
.262   X   10^ 
.022   X   10^ 
.093    X   10^ 
1.056    X   10^ 

1.05 

1.28 

960 

-512 

1.75 

1.06 

1325 

-147 

2.45 

1.02 

1785 

313 

3,50 

1.00 

2500 

1028 

7358 

1.906   X   10^ 

ge: 

w  =   1472 

• 

log 

sinh 

\ 

s 

1 

.2293 
.0982 
.4457 
.7597 
.2186 

Sum: 

Average 

^h 


log   r  -   log   sinh  b   s         (w-w)    x  prec. 


-.1796 
-.1440 
-.2443 
-.2386 
-.2766 

-1.0831 

=   -    .2166 

-127 


.4343   X   1.906   x   10 


123 

74 

36 

-    75 

-285 

-127 


=   -1.54    X    10 


-4 


Let      Ab   =   -1.5  X   IC 


-4 


Then 
-4, 


AA   =    .4343    (1.5    X    10      )(1472)    =    .0969 


The    corrected   values    of     b      and      /      are   b   =    .00125 

A   =   -.1206 
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s 

r  obs. 

r  calc. 

400 

.59 

.393 

753 

.90 

.820 

1250 

1.50 

1.73 

1750 

3.32 

3.34 

2500 

8.75 

8.60 

Using  log  a  =  A,  this  gives  a  =  .758.  The  final  formula  is 

r  =  .758  sinh  .00125s 

°^     '  r  =    .758  slnh  (   ^g^  ) . 

The  calculated  values  of  r  are  compared  v;lth  the  observed 
values  in  the  following  table: 

per  cent  error 

+  1.3 
-8.9 
+8.8 
+  .6 

-1.7 

The  graph  is  indicated  in  fig  2. 
Derivation  of  the  equations. 

To  minimize  (4)  set  its  partial  derivatives  with  respect  to 

A   and  b   equal  to  zero.   The  follov/ing  equations  are  obtained* 

(8)  Z_.  (log  r^  -  A  -  log  sinh  bs  )  =  0 

(9)  ^^^1  coth  bs .  (log  r.  -  A  -  log  sinh  bs.)=  0 
From  (8) 

(10)  A  =  i    /    (log  r   -  log  sinh  bs  ) 

n    ''  •■  ■  1  1 

In  (9)  the  "weights"  s.  coth  bs.  are  aprroxlmated  by  w.=s.    .  , 

1  X  xicot/n  D^  s ,  . 

O  1 

Solving   for      A      in    (9) 

wA  =  —     y     v/.(log  r.    -   log    sinh  bs.) 
n    ^        1        "      i  ^  1  ^ 

Eliminating  A   between  this  equation  rand  (10)  the  equation  for 
b   is  as  follows: 

(11)  /  (w.  -  ■.  )  (log  r.  -  log  sinh  bs.)  =  0 
Letting  b  =  b   +  A  b,  then 

(12)  log  sinh  bs.  =  log  sinh  b  s.  +  s.  coth  b  s, . /:^b(  .4343) 
approximately.  Substituting  (12)  in  (11)  and  solving  for  {S^h, 
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/    (w.    -   w)(log  r.    -   log   sinh  b^s    ) 
(13)  Ab  =  -^=-    ^  ^  °    - 


.4543«    2^  (w.    -   w)^ 
since   w.    =  s.    coth  b   s.    and   also    y    w.  (w.    -   w)    =   y     (w.    -  w)'^* 
Let   the   first   approximr:  tion   to      A      be 

Substitute    (12)    in    (10)    and    let     A      =  A      +^k,   we    find 

^A  =   -  -i  X  w.  /^-  b(.4343) 
n    — "      1 

or 

(14)  A  A   =   -.4343  w-Ab. 

Summary 

This  paper  demonstrates  how  to  fit  a  hyperbolic  sine 

curve  to  a  set  of  experimental  data.   A  least  squares  solution 

leads  to  a  pr.ir  of  transcendental  equations.  These  arc  solved 

by  iteration.   An  example  is  worked  out  which  shows  how  closoly 

this  law  fits  the  data  taken  from  the  Creep  Rate  Book. 
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